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Abstract. New Ha emission-line observations of the Virgo cluster galaxy NGC 4438 are presented. Fabry-Perot 
interferometry data at an effective angular resolution of ~ 2" are used to map the kinematics of the ionized gas in 
the galaxy. For the first time we obtain a velocity field covering a large area in NGC 4438, much larger than that 
deduced from previous Hi and CO maps. The kinematics of the extended, low surface brightness Ha filaments to 
the West of the galactic disk is discussed. We report on the discovery of a northern Ha structure which is clumpier 
than the other filaments. Evidence for multiple spectral components through the data-cube are presented in a 
nuclear shell and in the approaching half of the disk. The role of VCC 1040, a dwarf elliptical galaxy located to 
the South of NGC 4438, is presented to investigate the origin of a small-scale stellar tail of NGC 4438. It could 
be due to a minor tidal interaction between the two galaxies. 

Key words. Galaxies: individual: NGC 4438 - Galaxies: clusters: individual: Virgo - Galaxies: interactions - 
Galaxies : kinematics and dynamics - Galaxies: velocity field - Instrumentation: interferometric - Techniques: 
high angular resolution 


1. Introduction 

The role of environmental effects is crucial for the evo¬ 
lution of galaxies. Galaxies undergo several different pro¬ 
cesses that directly affect their interstellar medium, mod¬ 
ifying their morphology and perturbing their dynamics, 
particularly while crossing high density environments such 
as galaxy clusters. Well-known examples of such processes 
are tidal perturbations between galaxies (Moore et al. 
1996, Gnedin 2003) and ram pressure stripping exerted 
by the hot intracluster medium (hereafter ICM, Gunn & 
Gott 1972). 

The Virgo cluster galaxy NGC 4438 (VCC 1043) is 
a prototype of a disk interacting within a dense environ¬ 
ment. Its main properties are listed in Tab. ^ It is one 
of the closest objects to the cluster core at only ~ 1° in 
projection (or ~ 280 kpc) from M87 and ~ 28' (140 kpc) 
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from M86. Its morphology is the most perturbed among 
the Virgo galaxies, showing prominent stellar tails from 
North-East (NE) to South-West (SW) of the disk and a 
huge dust lane West of the disk (Fig. [TJ , molecular and 
neutral gas clouds displaced up to ~ 1' (4.7 kpc) West 
of the disk (Combes et al. 1988; Cayatte et al. 1990), 
very extended (up to ~ 2' or 9.4 kpc) ionized optical and 
X-ray emitting gas filaments (Kenney et al. 1995, here¬ 
after referred to K95; Machacek et al. 2004) and an elon¬ 
gated region of 1.4 GHz continuum emission (Hummel et 
al. 1983). It is one of the most Hi deficient disks in the 
cluster (Cayatte et al. 1990, 1994). Numerical simulations 
succeeded in explaining the NE and SW tidal tails forma¬ 
tion by a past high-speed tidal encounter with NGC 4435, 
a northern companion (Combes et al. 1988). While a role 
for ram pressure stripping by the ICM was proposed to 
explain some observations (Kotanyi et al. 1983, Cayatte 
et al. 1990, Keel & Wehrle 1993), K95 speculated that 
the off-plane filamentary components could be regarded as 
consequences of an ISM-ISM interaction with NGC 4435. 
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In addition to this strongly disturbed large-scale mor¬ 
phology, Ha+[Nn] WFPC2-H5T imagery (Kenney & Yale 
2002) and Chandra X-ray observations (Machacek et al. 
2004) have revealed gaseous nuclear bipolar shells. The 
origin of these nuclear shells could be due to an AGN 
(Machacek et al. 2004). 

Optical Fabry-Perot (FP) interferometry is quite ap¬ 
propriate to study the kinematics of NGC 4438 since the 
ionized gas is the only gaseous component covering both 
the inner disk regions as well as the outer ones. As part 
of a survey dedicated to the study of the kinematics of 
Virgo Cluster galaxies (Chemin et al. 2004, Chemin et 
al. in prep.), FP interferometry observations of the Ha 
emission-line in NGC 4438 are used to map its complete 
kinematics for the first time. Section [21 presents the high 
resolution velocity field of the galaxy, which reveals the 
kinematics in the disk, in the nuclear shells and along the 
external filaments. In Section 0 the kinematics of the fila¬ 
ments as well as the role of ram pressure stripping in their 
evolution are discussed. Section [S] also proposes that a mi¬ 
nor interaction with the dwarf elliptical galaxy companion 
VCC 1040 could be at the origin of morphological distur¬ 
bances observed in the southern region of NGC 4438. 

To be consistent with K95, a distance to the Virgo 
cluster of 16 Mpc is adopted throughout the article. 

2. Observation and results 

2.1. Data acquisition and reduction 

One set of observations of NGC 4438 was performed in 
April 2002 at the 3.6m European Southern Observatory 
telescope equipped with a focal reducer, a scanning FP 
interferometer and the photon counting camera FaNTOmM 1 
(Gach et al. 2002). Table [21 describes the characteristics of 
the observations. 

A FaNTOmM Ha emission-line observation obtained 
at the 1.6m Observatoire du mont Megantic (OmM) 
Telescope was also used to map a large field of view around 
NGC 4435 and NGC 4438 (see Tab. [21. The observation 
is used to search for emission-lines outside the disk in the 
large-scale tails of NGC 4438. Although the OmM data- 
cube has low signal-to-noise ratio, the detection of ionized 
gas was successful in the disk and the brightest parts of 
the filaments, as well as at three locations represented by 
white squares in Fig. [T| (see also Tab. El. Only the spectra 
at these three locations are presented for this observation 
(see section El- 

2.2. Ha Fabry-Perot emission-line and velocity maps 

The basic pre-processing data reduction steps of raw FP 
data-cubes in order to obtain a calibrated data-cube have 
been described elsewhere (see e.g. Chemin et al. 2003). 

1 FaNTOmM stands for Fabry-Perot de Nouvelle 

Technologie de l’Observatoire du mont Megantic. See 

http://www.astro.umontreal.ca/fantomm for more details. 


Table 1 . Basic properties of NGC 4438 (VCC 1043). 

Right ascension (J2000) 12" 27 m 45.6 s 

Declination (J2000) +13° 00' 32" 

Type 1 Sb (tides) 

Nucleus 2 LINER 1.9 

Distance 16 Mpc 

Linear scale 78 pc arcsec -1 

Heliocentric systemic velocity 3 v sys 71 km s _1 

Inclination 4 i 87° 

Disk position angle 3 P.A. 29° 


Morphological type from Binggeli et al. (1985) 
(2) Ho et al. (1997) 

^ v S ya and P.A. from Kenney et al. (1995) 

(4) i from LEDA 


Table 2. Observational setups. 


Telescope 

ESO 3.6m 

OmM 1.6m 

Observation date 

2002, Apr., 6 th 

2004, Feb., 24 th 

Field of view 

3.'6 

13.'7 

Data-cube size 

512 x 512 x 24 

512 x 512 x 48 

(X,Y) pixel size 

0^42 

l"61 

Interference order @ Ha 

793 

765 

Free spectral range @ Ha (km s — ) 

378 

392 

Channel width (km s~ ) 

15.75 

8.16 

Scanning wavelength (A) 

6564.3 

6561.7 

Interference filter central A (A) 

6567.7 

6565.0 

Interference filter FWHM (A) 

12 

30 

Exposure time (s) 

9360 

14400 

Average seeing 

00 

2'.'5 


Once calibrated, the integrated emission of the galaxy cor¬ 
responds to the sum of the signal of all channels above a 
continuum threshold. The threshold is chosen so that 15% 
of all the channels are under the continuum level, and the 
remaining 85% channels are consequently due to the Ha 
emission-line. We call this emission the total integrated 
emission. The pixel velocity corresponds to the barycen- 
tre of the emission-line. The case of a peculiar region for 
which some of its pixels exhibit two lines is described be¬ 
low. For these pixels, the velocity is not unique and a spe¬ 
cific treatment must be done (see “Southern disk-Filament 
F4 region” paragraph). 

Figure^) and FigureElpresent the total integrated Ha 
emission-line and the velocity field of NGC 4438 obtained 
using a data-cube appropriately smoothed to effective res¬ 
olutions between l'/26 and 2'.'52 (corresponding to 3 and 6 
pixels FWHM) in order to increase the signal-to-noise ra¬ 
tio in the field-of-view. Our Ha map is in good agreement 
with known deep Ha imagery (K95; Kenney & Yale 2002; 
Gavazzi et al. 2003). The Ha velocity field is in good agree¬ 
ment with the few long-slit measurements (K95) where 
both studies coincide. The most striking feature is the fil¬ 
amentary morphology extending West and South up to 2(1 
(or 9.8 kpc) from the nucleus and on which are superposed 
a few brighter clumps, probably Hii regions. Some of these 
structures have already been described (K95, but see also 
Keel & Werhle 1993) and similar names as in K95 for the 
filaments will be used hereafter, as shown in the bottom- 
left corner of Fig. Q]d. In the following, the kinematics of 
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Fig. 1. (a). A view of NGC 4435, NGC 4438 and VCC 1040. The optical map is a 13-band image taken from the 
GoldMine database (Gavazzi et al. 2003). Contours display the approximate coverage of the total Her emission. White 
squares correspond to the positions of selected spectra from the Megantic FP observation (see 1 12.21 Tab. Inland Fig. El)■ 
The coordinates are in J2000. Notice the large dust layer on which the Ha emission superposes to the North and West 
of the disk of NGC 4438. A dashed box corresponds to the field of view seen in the right panel, (b). Contours and 
grey-scale (displayed using a logarithmic stretch) map of the total integrated Ha emission. The first contour is at a 
higher intensity than the minimum level of the grey-scale image. The bottom-left, reduced-scale cartoon indicates the 
nomenclature used for the filaments. 
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Fig. 3. Ha Position-velocity diagrams through different 
regions of the continuum-free data-cube. Top : along the 
major axis of the disk (P.A. = 29°). Bottom : along the 
Filament F4 (P.A. = 180°). A logarithmic grey-scale and 
contoured stretch is used for the Ha emission-line. A red 
solid line corresponds to the velocity field values (Fig. [3). 
The centre and orientation of the slices are listed in Tab. [3 
The effective angular resolution of the data-cube used for 
the diagrams is 0"84 (or 2 pixels FWHM). 



Fig. 4. Integrated emission (a) and velocities (b) of the 
most redshifted line for peculiar pixels in the dashed box 
region of Fig. |3 (see section l?~2l paragraph “Southern disk- 
Filament F4 region” for details), (a). The grey-scale pixels 
display the emission integrated under the most redshifted 
emission line of these peculiar pixels and contours display 
the total integrated emission of the galaxy (already seen 
in Fig. Et)- (b). Same velocity field as in Fig. [2] but this 
time showing the red-shifted velocity component for the 
peculiar pixels. Notice that these redshifted lines have ve¬ 
locities in contrast with that of the approaching half of 
the disk but in agreement with those of structures outside 
the disk and Filament F4 (encircled by dashed lines). 


the different components of the galaxy are described: disk, 
nuclear shells, filaments and new found structures. 


Disk. In the disk, the gas emission is dominated by 
few Hil regions South and North from the nucleus and by 
a nuclear shell (see next paragraph). No attempt to derive 
an optical rotation curve was undertaken due to a highly 
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Table 3. (Top) Position-Velocity slice parameters along 
the major axis and Filament 4 from the ESO FP- 
observation. (Bottom) Position of selected spectra in the 
OmM FP-observation. 


ESO : 



Slice 

Centre coordinates 

P.A. 


(J2000) 

n 

Major axis 

12 h 27 m 45=698 ±13° 00' 31738 

29 

Filament 4 

12 h 27 m 44=997 ±12° 59' 497 57 

180 

OmM : 
Spectra 

Centre coordinates 

Velocity 


(J2000) 

(km s _1 ) 

South (a) 

12 h 27 m 38!887 ±12° 58' 327 73 

-35 

North (b) 

12 h 27 m 46!410 ±13° 02' 29'.'27 

-60 

North (c) 

12'* 27 m 48!000 ±13° 02' 45"00 

-200 


inclined disk, obscured by dust lanes, that would cause 
large uncertainties in the circular velocities. The average 
minimum velocities are down to —150 km s -1 in the south¬ 
ern approaching half of the disk. Our results are consistent 
with the CO observations which show an inner molecular 
disk with velocities lower than —100 km s -1 (see Figs. 
10 and 12 of K95, see also discussion in K95 about the 
Combes et al. observations). 

Figure0represents the position-velocity diagram along 
the major axis of the galaxy (P. A. = 29°) and shows a clas¬ 
sical rotating disk. The kpvslice task of the KARMA pack¬ 
age (Gooch 1995) was used to make the position-velocity 
diagram. The pseudo-slit width is between 3 and 4 pixels 
(i.e. between 173 and l'/7) with a pixel weighting inversely 
proportional to the distance of the slice. For clarity rea¬ 
sons, a red solid line, which represents the velocity field 
values (Fig. [3), is drawn to follow the main interference 
order since more than one order is shown here. Lines are 
broader in the centre than in any other regions of the ma¬ 
jor axis. 

Nuclear shells. Kenney & Yale (2002) found ev¬ 
idence for a bipolar nuclear shell in NGC 4438 from 
Ha+[Nil] WFPC2-PST images. The NW and SE shells 
are intrinsically asymmetric in luminosity and length, in 
both optical and X-ray emission-lines and radio contin¬ 
uum (Kenney & Yale 2002, Machacek et al. 2004). Using 
an emission-line velocity dispersion measurement along 
the western edge of the NW shell, Kenney & Yale (2002) 
claimed that the shells are caused by a nuclear outflow in 
NGC 4438. Figure ^ shows that the NW shell is one of 
the brightest structure that is observed in the disk (R.A.~ 
12^ 27 m 45.7 s , Dec.~ +13° 00' 33") while the SE shell 
is indeed fainter (R.A.~ 12^ 27 m 46.5 s , Dec.~ +13° 00' 
26"). Figure [3 shows that most of the gas motions in the 
SE shell is at u 0 bs ~ 75 ± 5 km s -1 , which is close to the 
systemic velocity (71 km s -1 ). At the tip of the shell, the 
velocities reach ~ 115 km s^ 1 . 

In Figure^ we focus on the NW shell kinematics. The 
shell has radial velocities increasing from ~ —60 km s -1 to 
~ 140 km s -1 (from SW to NE of the shell, respectively). 
We use here the high resolution Ha + [Nil] WFPC2-P5T 


images (Kenney & Yale 2002) taken from the HST archive 
and reduced as explained in m The velocities are close 
to the systemic velocity, principally along the apparent 
semi-major axis of the shell. As seen in spectra extracted 
from the full angular resolution data-cube and selected at 
different positions along the shell (right-panel diagrams of 
Fig. 0), the profiles are relatively wide in the shell. The 
restricted free spectral range of 380 km s -1 of the ESO- 
FP observation does not allow us to derive the velocity 
dispersions because the bases of the profiles are missing. 
The spectra show several velocity components. The major 
component in each spectrum is at a velocity given by the 
velocity field, close to the systemic velocity of the galaxy 
of ~ 70 km s -1 . A redshifted (blue-shifted, respectively) 
component is detected at velocities up to ~ 150 km s -1 
(down to ~ 50 km s” 1 ) relative to the systemic velocity. 
The blue-shifted component is detected in diagrams 1 to 3, 
but hardly in diagram 4. The reader should be aware that 
these additional components are known ±380 km s -1 , i.e. 
the ESO interfringe, because there is no direct way to de¬ 
termine the interference order at which they are emitting. 
Determining their true velocity would need larger spec¬ 
tral band observations than ours. Line splittings in nu¬ 
clear shells often indicate an expanding bubble (Veilleux 
et al. 1994). This is here the first direct evidence that 
the shell is expanding. An accurate determination of the 
intrinsic outflow velocity and velocity dispersion would re¬ 
quire a coupled kinematical and geometrical model of the 
shell (see e.g. Veilleux et al. 1994, for NGC 3079), which 
is beyond the scope of the article. 

Filament FO. This northern filament was not previ¬ 
ously described in K95 although present in their data. It is 
located at 65" (5.1 kpc projected) North from the nucleus 
and is mostly traced by a string of Ha knots (probably Hii 
regions) embedded in faint diffuse gas, running nearly par¬ 
allel to the EW direction. Its morphology thus differs from 
that of the other filaments, which are mostly only diffuse. 
It has a length of 45" (3.5 kpc) although it could be longer 
toward the SW and connected to Filament FI. Compared 
with a B-band image of the galaxy (Fig. QJ, F0 exactly 
superposes on the northernmost part of the large western 
dust layer. Its velocities (~ 80 km s -1 ) are of the same 
order as the systemic velocity of NGC 4438 but admit a 
maximum of ~ 110 km s _1 at its western tip - the Hii 
clump at R.A.=12' 1 27 m 44.4 s , Dec.=±13° 01' 36". The 
velocity of this latter clump is consistent with that of a Hi 
counterpart at 120 km s -1 (Cayatte et al. 1990) as well 
as the one obtained with the long-slit observation (K95, 
their P.A. = 0° slit). A comparison of our Ha map with 
the Chandra X-ray image (Machacek et al. 2004) shows 
that F0 has no X-ray emission counterpart. Its evolution 
appears slightly different from that of the other filaments. 

Filament FI. A velocity gradient is observed along 
part of the filament over a distance of 20" (1.6 kpc) from 
250 km s -1 down to about the systemic velocity of 71 km 
s -1 (from East to West respectively). One may notice that 
the western tip of FI actually presents two regions with 
distinct velocities: a northern part at u Q bs ~ 130 km s -1 
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(R.A.~ I2 h 27 m 41.6 s , Dec.~ +13° 01' 10") and a south¬ 
ern one at the systemic velocity. In fact, this northern 
part has a velocity comparable with that of the western 
extremity of F0 (within 20 km s -1 ) and could be a genuine 
extension of F0. Combes et al. (1988) reported wide CO- 
line profiles at a velocity of ~ 85 km s^ 1 in the molecular 
cloud located at V (4.7 kpc) NW from the nucleus. Taking 
into account the resolution of 23" for the CO observations, 
it is probable that these two regions at significantly differ¬ 
ent velocities could not be resolved by the IRAM beam 
and thus artificially enlarged the CO emission-lines. These 
latter are then nearly centered at the velocity of the high¬ 
est surface brightness component. 

Filament F2 - This filament also presents a velocity 
gradient, with velocities between ~ —15 km s -1 at its 
western tip and ~ 200 km s _1 at the apparent connection 
with the disk. Compared with the largest neutral hydrogen 
concentration detected in the VLA observations (Cayatte 
et al. 1990), the Ha and Hi velocities are very similar for 
the regions of overlap (60-80 km s -1 ). 

Filament F3.~ It is composed of three distinct, nearly 
parallel substructures F3a, F3b and F3c. F3b represents 
the main body of F3 and exhibits an average velocity of 
Fobs ~ 20 km s -1 occurring in its brightest parts, although 
lower velocities down to ~ —30 km s -1 are also observed 
along it, as well as higher ones (up to ~ 70 km s -1 ) at its 
eastern tip, close to F4. F3b actually presents wide profiles 
(up to 165 km s -1 FWHM, corrected for the instrumental 
Airy profile) which could be due to turbulent motions. 

F3a has negative velocities (~ —25 km s -1 ) to its NE 
side and null velocities (on average) on the other side. A 
new, short filament (F3c) is detected (~ 15 ", 1.2 kpc), 
barely seen and not discussed in K95. It appears to be 
linked to F3b by a ~ 5" (400 pc) ionized gas structure. 
F3c appears to be an extension of a bright warped stellar 
structure (see Fig. EK) and mostly shows velocities of —20 
km s , but with a drop down to —80 km s _1 . 

Filament Ff. - F4 appears to be the longest filament 
with a projected length of at least 1.'5 (7 kpc). Its first 
part originates in the southernmost side of the disk and is 
along the NS direction, running parallel but clearly offset 
from a stellar tidal tail (see - I!..'ll and Fig. 0i). Its second 
part is directed along the NE-SW direction and is almost 
parallel to F3. 

F4 is the most blue-shifted filament w.r.t. the systemic 
velocity, presenting an average and nearly constant veloc¬ 
ity of Fobs ~ —85 km s^ 1 . F4 may contribute to part of 
the kinematics in the southernmost parts of the disk but 
it is obviously less blue-shifted than further inside in the 
disk, where velocities down to —150 km s -1 are observed. 
It also exhibits a velocity drop (down to ~ —135 km s^ 1 ) 
occurring in its bending region. Ionised gas at Fobs ~ —85 
km s -1 is also observed close to the bending region but 
outside F4. The clump at the tip of F4 is at ~ —5 km s -1 
and therefore appears kinematically independent from the 
whole filament. 

Southern disk-Filament FJ f region. - The dashed 
box of Fig. 0 corresponds to a particular region, where F4 


appears to meet the SW half of the disk. Some of its pix¬ 
els exhibit two emission lines (Fig. 0). For those pixels, a 
barycentre computed over the whole profile, as done for all 
other pixels of the field-of-view that have only one spectral 
line, would give incorrect velocities due to the presence of 
the two competing lines. Two velocities need to be com¬ 
puted for them. A blue-shifted component has a velocity 
comparable with that of neighbouring pixels that show 
only one spectral component (< —100 km s -1 ). It allows 
to follow a kinematic continuity in the SW approaching 
half of the disk. Its pixels are thus shown on the velocity 
map of Fig. 0 Figure 0 shows the pixels of the second, 
redshifted component (up to ~ 100 km s -1 ). A disconti¬ 
nuity is clearly seen with respect to the surrounding pixels 
showing only one emission line. Its kinematics is hard to 
reconcile with the rotating, approaching, blue-shifted pat¬ 
tern of the SW half of the disk and of Filament F4. It 
does not suggest that the emission comes from the disk 
or Filament F4, unlike the blue-shifted line. Its velocity 
is in better agreement with that of diffuse Ha structures 
outside the disk and F4 (see circles in Fig. QJ. The pix¬ 
els seem to connect these latter different redshifted struc¬ 
tures, as if all the redshifted gas would belong to a same 
morphological structure. This latter could be an exten¬ 
sion of Filament F2 to the East. The amplitude of the 
most redshifted line tends to be brighter than the blue- 
shifted one. However the transmission curve of the narrow 
interference filter may bias this claim (the filter transmis¬ 
sion being higher in this part of the spectral range than 
in its blue part) so that it is hard to discriminate which 
component really dominates the total integrated emission. 
The integrated intensity map of the redshifted line is rep¬ 
resented in grey-scale in Fig.0 with contours showing the 
total Ha intensity. This local map should only be used to 
locate the peculiar pixels. Further observations would be 
needed to confirm our new finding and to accurately map 
the intensity of the most redshifted component. 

Evidence for new structures. An isolated clump 
is detected West of F0 (R.A.=12 ,t 27 m 42.5 s , Dec.=+13° 
01' 44") with a radial velocity of 145 ± 5 km s _1 . It is 
barely seen in the Ha image of Kenney & Yale (2002) 
and no other molecular gas or Hi counterparts have been 
observed at this location. The Ha image of Kenney & Yale 
(2002) also shows two other isolated knots located at 60" 
and 65" NE from the nucleus and close to the eastern 
edge of their field of view. These knots do not present 
Ha emission in the FP data-cube but typical continuum 
profiles and thus they could be foreground stars or more 
distant ojects. 

To the West of the field of view, in a region centered 
at R.A. ~ 12 h 27 m 39 s , Dec. ~ +13° 00' 00", a very faint 
and patchy emission is observed at velocities between —80 
and 20 km s -1 . It is possible that sky emission-line resid¬ 
uals contaminate this extremely low surface brightness re¬ 
gion. Nevertheless, most of this emission can be seen in 
the Kenney & Yale (2002) Ha image, although not dis¬ 
cussed, as well as in the Ha image of GoldMine (Gavazzi 
et al. 2003). Therefore the filamentary structure appears 
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much more complex and extended towards the West than 
the previously discussed linear filaments. 

To the East of F4, a new faint Ha structure is ob¬ 
served, with clumps at velocities between 70-90 km s -1 
(see e.g. the clump at R.A. ~ 12 h 27 m 46^053, Dec. ~ 
+12° 59' 54'.'06). These Ha clumps superpose on a short 
X-ray filamentary-like structure in the Chandra observa¬ 
tions of Machacek et al. (2004), although not discussed by 
them. The Ha clumps are separated from F4 by a dust 
lane (see Fig. Q. We think that we do not detect a fila¬ 
ment like the X-ray one because our FP profiles are at too 
low signal-to-noise in this region. For the present time, the 
origin of the Ha clumps and of the short X-ray structure 
is unknown : they could have a relation with F4, and/or 
with the nuclear outflow, or perhaps with the pixels that 
show a redshifted spectral line in the southern disk region 
(see previous paragraph). 

Three clumps outside the ESO-FP held-of-view have 
been detected with the OmM observations. Their spectra 
is shown in Fig. El and their velocities listed in Tab. 0 
Their positions are marked by white squares in Fig. ^ 
The lines observed here cannot be associated to night sky 
emission-line residuals because these latter always present 
very sharp wings and are relatively narrower than ionized 
gas emission-lines, which is not observed. The left panel 
(Fig. Eli) shows that ionized gas is detected in the SW 
stellar tail, in the continuity of F3b and with a line at a 
velocity of u 0 b s ~ —35 km s . This is comparable with the 
velocities observed along F3b. On broad-band images, the 
knot associated with the middle panel spectrum (Fig. EJj) 
appears as a system having two objects. The right-hand 
spectrum (Fig. Et) corresponds to a blue clump on the 
RGB image of the GoldMine database. These two latter 
clumps could be Hll regions associated with the northern 
stellar tail. Their velocities appear extremely blue-shifted, 
in contrast with the kinematics observed in the northern, 
receding half of the disk. However, because the velocity of 
any isolated clump is known ±390 km s^ 1 (i.e. the OmM 
observation free spectral range), it is reasonable to deduce 
true velocities of ±330 km s -1 and ±190 km s -1 for them, 
when adding an interfringe to the measured velocities. In 
this case, their velocities would be in better agreement 
with the kinematics in the North. Their velocities need to 
be confirmed by other spectroscopic measurements. 
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Fig. 6. Fabry-Perot continuum-free spectra of selected 
regions (see coordinates in Tab. El and white squares in 
Fig. □. The spectra are averaged in a box of 8" x 8" 
around centres of Tab. El The flux unit is arbitrary. The 
dashed lines indicate the systemic velocity of the galaxy. 


3. Discussion 

3.1. Global kinematics of the filaments 

Evidence for a high-speed tidal interaction (~ 100 Myrs 
ago) with the northern companion NGC 4435 at the origin 
of the prominent NE and SW stellar tails and of part of the 
molecular gas above the galactic disk were first presented 
in Combes et al. (1988). Furthermore, K95 proposed that 
the off-plane gaseous distribution could be due to an ISM- 
ISM collision with NGC 4435 and that the gas should be 
re-accreted into NGC 4438 after the collision. 

Numerical simulations of galaxy tidal interactions or 
mergers (e.g. Hibbard & Mihos 1995, Barnes & Hernquist 
1996, Mihos 2001) have shown that most of the mate¬ 
rial in stellar and gaseous tails remain bound to a parent 
disk, hence falling back towards it. There are evidence 
in the velocity held that gas in the filaments is dynami¬ 
cally bound to NGC 4438. Indeed, this study confirms that 
no extremely high velocity departures with respect to the 
systemic velocity are observed in the filaments (see also 
Combes et al. 1988, K95). The highest velocity amplitudes 
(in absolute values) are found in the disk. Moreover, the 
filaments seem to participate to the rotation of the galaxy, 
with redshifted (blue-shifted) gas to the North (South re¬ 
spectively) and with velocities close to the systemic value 
in parts of FI and F2 that roughly lie along the minor 
axis. Though F0 is obviously less redshifted than part of 
FI while it is the northernmost filament, this trend exactly 
follows the disk rotational pattern. The diffuse filaments 
of NGC 4438 are thus apparently bound to it and it is 
likely that they will be reaccreted towards the disk. 


3.2. Ram pressure stripping 

Signatures of ram pressure stripping (hereafter RPS) have 
already likely been observed in nearby groups and clusters 
galaxies in the form of linear structures, filaments, off- 
plane ring-like and comet-like morphologies of either ion¬ 
ized or neutral gas or radio continuum (Gavazzi et al. 2001, 
Vollmer et al. 2000, Bureau & Carignan 2002, Vollmer et 
al. 2004a). In the Virgo cluster, there exist observational 
indications of ongoing or recent RPS. The most probable 
cases are among the most Hi-deficient galaxies, like NGC 
4388 (Veilleux et al. 1999, Yoshida et al. 2002, Vollmer 
& Huchtmeier 2003) or NGC 4522 (Vollmer et al. 2000, 
Kenney et al. 2004, Vollmer et al. 2004a). A case where 
a RPS could happen with a gravitational interaction has 
also been presented for NGC 4654 (Vollmer 2003). 

It has been claimed that a tidal interaction and an 
ISM-ISM collision with NGC 4435 were probably respon¬ 
sible for the displacement of Hi and part of the CO gas 
outside the galaxy plane of NGC 4438 and for the for¬ 
mation of the ionized gas filaments (Combes et al. 1988, 
K95). NGC 4438 is yet a good candidate for an active 
RPS by the Virgo ICM (Kotanyi et al. 1983, Cayatte et 
al. 1990, 1994, K95, Vollmer et al. 2004b). The fact that 
F2, F3, parts of F0 and FI and the second half of F4 are 
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all directed along the NE-SW direction could be an indi¬ 
cation that a common stripping process is responsible for 
their shaping and high degree of collinration. 

A key-point for a RPS hypothesis is to explain why no 
gas is observed at high velocities relative to the galaxy in 
the Ha filaments. This observational fact is not expected 
because of the rapid 1-o-s relative motion of NGC 4438 
with respect to the cluster (~ 1000 km s _1 ). If high veloc¬ 
ity gas really exists, which reasons could explain its non¬ 
detection ? Firstly, we recall that our restricted spectral 
ranges prevent from detecting any extremely high velocity 
departures w.r.t. the systemic velocity. If present, detect¬ 
ing it would need other observations performed at other 
wavelengths or with a larger spectral band. However, the 
Fabry-Perot spatial distribution of Ha matches well the 
one obtained from deep Ha + [Nil] images, showing that 
no gas is missed by our observation. In other words, it 
seems implausible that a large amount of gas at high 1-o-s 
velocity exists in the optical filaments. Then, numerical 
simulations of RPS could help to answer this question. 
Major improvements on the modelling of RPS have been 
done recently (Abadi et al. 1999, Quilis et al. 2000, Schulz 
& Struck 2001, Vollmer et al. 2001) and models succeeded 
in reproducing kinematical properties and morphological 
anomalies of cluster galaxies (e.g. Vollmer et al. 2004b). 
For instance, the case of NGC 4388 is particularly in¬ 
teresting because it is another highly inclined galaxy in 
the cluster core. A RPS model for NGC 4388 (Vollmer 
& Huchtmeier 2003) has shown that most of the gas re¬ 
mains bound to the galaxy and that gas particles with the 
highest velocities w.r.t. the galaxy appear at high distance 
from the disk (few tens of kpc). The important parame¬ 
ters of such simulations are the ICM density, the galaxy 
velocity relative to the ICM (2000 km s -1 for NGC 4388) 
and the inclination angle of the galaxy disk to its orbital 
plane. Although NGC 4438 is orbiting in the Virgo clus¬ 
ter under different conditions than NGC 4388 and could 
have been considerably disturbed by a massive compan¬ 
ion, one could expect that most of its gas resides close the 
galaxy systemic velocity within the first tens of kpc from 
the centre, as is observed. If a small fraction of gas with 
high relative velocity exists, it could be now located fur¬ 
ther away from the nucleus, outside the field of view. A 
good test would be to search neutral gas in a very large 
field-of-view around the galaxy. 

Perhaps another good test for an ongoing RPS scenario 
would be to study the stellar kinematics of the southern 
region, where filaments F3 and F4 appear to superpose 
on large and small (respectively) stellar tails. Since RPS 
only acts on ionized and neutral gas and not on the stars, 
it is expected that a kinematical decoupling has devel¬ 
oped between gas and stars. It should however be remem¬ 
bered here that simulations of galaxy collisions have al¬ 
ready shown velocity differences reaching ~ 50 km s -1 
between stars and gas in tidal tails (Mihos 2001), though 
their global evolution remains coupled (re-accretion on the 
parent galaxy disk, ...). We do not know the amplitude 
of an hypothetical velocity difference between stars and 


gas in the concerned regions of NGC 4438 that would 
be caused by a tidal interaction alone. Numerical mod¬ 
els of RPS coupled with a tidal interaction should help to 
disentangle their respective role on the evolution of the 
filaments. 

3.3. A perturbed region to the South 
3.3.1. The peculiar Filament F4 

Filament F4 is probably the most intriguing observed fila¬ 
ment. A R-band image taken from the GoldMine database 
and F450W, F814W WFPC2 -Hubble Space Telescope im¬ 
ages taken from the HST archive are used to illustrate 
the tight relation that exists between stars, gas and dust 
around Filament F4 and the outer southern part of the 
disk (Figs. 0 to El). The WFPC2-R5T images were re¬ 
duced using the STSDAS package in iraf. Warm pixels 
and cosmic ray hits were removed using the warmpix and 
crrej tasks and the final mosaics were created using the 
wmosaic task. Cosmic rays were also treated using the 
lacos-im procedure (van Dokkurn 2001). We also built a 
F450W-F814W (~ B — I) colour index map from the ra¬ 
tio of F450W to F814W (Fig. 0j) in order to follow the 
dust lanes and star formation regions across and outside 
the galactic disk. Basically, redder colours appear here as 
darker shades and correspond to dusty regions (dust lane 
associated with the disk and filaments, and the western 
dust layer) while bluer colours appear as lighter shades 
and can correspond to star formation sites. The colour in¬ 
dex image principally shows that the dust is along a promi¬ 
nent ring-like structure in the inner parts of the galaxy 
and also present outside the disk to the West and South. 
Hu regions along the filaments have colours as blue as 
WF450W — m.F8i4W = —1-3 mag. 

We notice here the presence of an irregular-shaped 
object with a colour of up to 1.2 mag, lying along the 
edges of the prominent dust lane and Filament F3b (R.A. 
~ 12 h 27 m 42!3, Dec. ~ +12° 59' 35"2). It is nearly 
aligned with the bright warped stellar structure of the 
southern edge of the disk and Filament F3c (Fig. 0i). 
The only referred object closest to this location is called 
2MASXU1227422+125935 and is classified as a galaxy of 
diameter 1.'2 in the 2MASS database (Jarrett et al. 2000). 
We deduce an upper limit “diameter” of ~ 6" from the 
F814W-HST image, which is considerably lower than the 
referred value. No gas is detected within this object and 
because its redshift is unknown, it is not yet possible to 
determine whether it is related to NGC 4438 or whether 
it is a background galaxy. 

As seen in Figs. 0i and b, a first half of F4 lies along 
the eastern edge of a stellar tail (hereafter Tl) and exactly 
superposes on a dust lane. Then F4 bends towards the SE 
from the tip of Tl. The northernmost part of Tl seems 
to arise from the inner regions of the disk. A lower limit 
for the length of Tl is V (4.7 kpc). It is shorter than the 
two large-scale outer tails, which extend on distances of at 
least 4.'1 (or 19 kpc) towards the NE and 3.'1 (or 14.5 kpc) 
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Ellipse semi major-axis [arcsec 1 / 4 ] 

Fig. 8. Isophotal analysis of VCC 1040 from a F814W 
WFPC2-HST image. From Top to Bottom: intensity, el- 
lipticity, position angle (in degree), central X and Y co¬ 
ordinates (in pixel, full dots and open triangles respec¬ 
tively) variations as function of ellipse semi-major axis (in 
arcsec 0 25 ). The arrows indicate the radii of the isophotes 
displayed in Fig. M 


towards the SW. T1 is obviously crossed by a brighter and 
warped stellar tail (R.A.~ 12 h 27 m 44.5 s , Dec.~ +12° 59' 
55", Fig. 0i). Figure CJi shows that this bright warped 
stellar tail arises from the SE near side of the galaxy disk. 
T1 then lies straight towards the South, pointing in the 
direction of VCC 1040, a dwarf elliptical galaxy. It has 
a lower surface brightness than the stellar disk and the 
bulge/halo of NGC 4438. The tips of T1 and the brighter 
warped tail have blue colours (—0.2 mag. on average). T1 
is reminiscent of tidal stellar streams like those observed 
around massive spirals such as the Milky Way (Ibata et 
al. 2001a) or Messier 31 (Ibata et al. 2001b). The next 
section presents arguments to explain its formation by the 
presence of VCC 1040. 

3.3.2. A role for VCC 1040 ? 

VCC 1040 is a dwarf-nucleated elliptical galaxy and is the 
closest companion of NGC 4438. It is projected at only E5 
(7 kpc) from the nucleus of NGC 4438 (to be compared 
with 4.'5 or 21 kpc between NGC 4435 and NGC 4438) 
and ~ 25" (2 kpc) from F4. It has a redshift of —143 ± 
63 km s -1 (Conselice et al. 2001) which coincides within 
the errors with the velocity of the region where F4 bends 
and with the mean radial velocity of F4. Though dwarf 
elliptical galaxies are the most common type of galaxies 


found in the local galaxy clusters like Virgo (Ferguson & 
Binggeli 1994) and the probability of finding foreground or 
background dEs in the direction of the cluster core is high, 
their proximity in sky position and redshift likely makes 
VCC 1040 a genuine companion of NGC 4438. Its colour 
index B — R is 1.6 (K95) and indicates an old object (i.e. 
with an age >> 100 Myrs), which is usual for dEs (Rakos 
et al. 2001, Conselice et al. 2003, Jerjen et al. 2004), and 
no ionized gas is found in it. 

K95 have noticed the presence of several apparent stel¬ 
lar clump debris including VCC 1040, which could have 
been formed during the interaction between NGC 4435 
and NGC 4438. Tidal dwarf galaxy candidates are usually 
blue, gas rich, young and irregular-typed structures always 
observed along stellar tidal tails (e.g. Due et al. 1997), 
which properties do not match with those of VCC 1040. 
Morever the time-scale invoked to form a dE is several 
hundred Myrs (Kroupa 1998), which is shorter than the 
elapsed time since the interaction between NGC 4435 and 
NGC 4438 (~ 100 Myrs according to the Combes et al.’ 
simulations). It seems unlikely for VCC 1040 to be a tidal 
dwarf galaxy. 

Several elements show that the environment around 
VCC 1040 and VCC 1040 itself are perturbed. Firstly, an 
Ha emission is found between Filament F4 and the dwarf 
galaxy (Fig. |J|i). This emission actually appears to be an 
unusual ionized gas loop to the East of the F4 bending 
region. 

Secondly, the isophotal analysis of the dE light dis¬ 
tribution reveals its perturbed morphology, as shown in 
Figs. IH1 and El Ellipses of variable geometrical parameters 
(centre coordinates X0 and Y0, ellipticity and position an¬ 
gle of the isophotes) were fitted to a F814W WFPC2-H5T' 
image. A lower limit of the VCC 1040 radius is thus ~ 10" 
(780 pc), i.e. for the last accurately measured isophote. 
K95 found an outer radius of 14". At higher radius than 
10", the light distribution is partly contaminated by the 
diffuse stellar light of NGC 4438 and by dust to the N and 
NW outer parts of the dE, hence results are not shown for 
these radii due to too large uncertainties on the parame¬ 
ters. The graphs show that there is a substantial offset of 
the position of the photometrical centre as the radius in¬ 
creases and that the Y-offset of the outer ellipses is larger 
than the X-offset. As an example, the isophote with semi¬ 
major axis r ~ 8" 4 (last arrow) has its centre offset by 
~ 1"8 (140 pc) to the N-NW direction with respect to 
the innermost isophotes, within a ratio of xl°ffsct = 
Moreover, the galaxy is less rounder in its outer parts 
(e —> 0.2) than in its inner ones (e ~ 0.05 for 0"4 < r < 
1"), though passing through a significant peak at r=l"l (e 
= 0.4). Twists of the isophotes are also clearly detected. 
The P.A. indeed drastically varies between ~ 10° and 70° 
(for r > 0.5"). At a radius of ~ 0"25, the intensity bump, 
the ellipticity peak and the drastic change of photometric 
P.A. are likely due to the nucleus of the dwarf. Figure El 
displays some isophotal ellipses (which semi-major axis 
positions are flagged by arrows in Fig. El and illustrates 
well the measured asymmetry of the dE. 
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Fig. 9. Enlargement of an F814W WFPC2-HST image 
around VCC 1040. The arrows indicate different clumps 
surrounding the dwarf galaxy while ellipses indicate fitted 
isophotes of VCC 1040 to the image at the radii flagged 
by arrows in Fig. 0 The bottom-left insert is a focus of 
the white dashed box for another intensity scale to em¬ 
phasize on brighter clumps between the dE and the tip of 
the stellar tail Tl. 


Finally, numerous clumps are observed in the close 
vicinity of VCC 1040 as well as between the dwarf galaxy 
and the tip of Tl (Fig. 0. Their colour range is large, 
from blue to red clumps (—1 to 2 mag) indicating that 
they could be gas rich objects to old globular clusters. 
No gas was detected in them and their radial velocities 
are unknown so that their possible association with the 
southern region remains to be confirmed. Other knotty 
structures are found within the NE and SW large-scale 
stellar tidal arms (K95). It is surprising to find stellar 
clumps all around the dE, because they are not directly 
associated to the large-scale tails. Potential reservoir of 
stars for them could thus be VCC 1040 or Tl. 

Which mechanism could explain all these perturba¬ 
tions ? We emphasize here on a minor tidal interaction 
event that could have occurred (or perhaps that would 
be occuring) between VCC 1040 and NGC 4438, due to 
the presence of the dust lane and of the stellar tail (Tl) 
both pointing towards VCC 1040 and to the fact that the 
VCC 1040 redshift is comparable with velocities observed 
along the gaseous filament F4, and thus perhaps with the 
velocity of the tail. A comparison of the VCC 1040 tidal 

radius r t = R 3 ~ 7" (0.5 kpc) with the VCC 1040 
radius of r <; 10" (0.78 kpc) implies that a tidal strip¬ 
ping of its external parts by the gravitational forces of 
NGC 4438 could be ongoing. We have used here R ~ 7 kpc 
as the distance between the two galaxies, m oc 10 s M 0 as 
the mass of VCC 1040, which is typical of dwarf ellipticals, 


M oc 10 11 M 0 for NGC 4438 (Combes et al. 1988) 2 . Early 
simulations of the encounter with NGC 4435 (Combes et 
al. 1988) did not predict Tl but only the large-scale stellar 
asymmetry of NGC 4438 and one has to imagine another 
triggering event for its formation. The stream Tl and 
the stellar concentrations around the dwarf galaxy can be 
consequences of such a minor tidal event. Photometrical 
anomalies are often observed in dwarf early-type galaxies 
of the Virgo cluster (Binggeli et al. 2000, Barazza et al. 
2003) and it is still unclear whether they are due to in¬ 
trinsic or external perturbations to the dwarf galaxies. We 
think that the asymmetries measured in VCC 1040 are 
of the same nature as those measured in other dEs. An 
interaction with NGC 4438 could also explain the photo- 
metrical disturbance of VCC 1040, by stretching the outer 
parts towards NGC 4438, as measured from the isophotal 
analysis. 

This would not be the first time that a minor inter¬ 
action opposing two galaxies, one large and one dwarf, is 
observed in the Virgo cluster. Other possible cases have 
indeed been presented with the couples M49-VCC 1249 
(Sancisi et al. 1987, McNamara et al. 1994), NGC 4694- 
VCC 2062 (van Driel & van Woerden 1989, Hoffman et 
al. 1993) and M86-VCC 882 (Elmegreen et al. 2000). A 
detailed analysis of the stellar populations in the tails and 
in the galaxies would be very useful to test the minor in¬ 
teraction hypothesis. 

4. Conclusions 

New optical observations in the Ha emission-line of the 
Virgo cluster galaxy NGC 4438 were presented, as part of 
a large survey dedicated to the kinematics of Virgo spiral 
galaxies. Using Fabry-Perot interferometry data at an ef¬ 
fective angular resolution of ~ 2", the optical velocity field 
of the galaxy was mapped for the first time. NGC 4438 is 
one of the most disturbed Virgo spirals which has probably 
undergone a tidal interaction with the massive companion 
NGC 4435, and which is probably undergoing the effects 
of ram pressure stripping by the hot ICM in the cluster 
core. The Ha velocity field reveals the complex gas kine¬ 
matics within the galactic disk, the nuclear shells and the 
external filaments. Our conclusions can be summarized as 
follows: 

(1) The nuclear shell of NGC 4438 is likely expand¬ 
ing, as revealed by the multiple spectral components de¬ 
tected in the FP data-cube along the NW shell. Further 
observations of the nuclear region at higher spectral reso¬ 
lution and larger spectral coverage than our observations 
are needed to investigate the internal motions of the nu¬ 
clear outflow. 

(2) A new filament is observed to the North of the disk. 
It is short, clumpy and has no X-ray counterparts while 
the other filaments are much more extended, diffuse and 

2 The estimated value of the tidal radius is a lower limit 
because of the use of the projected distance of the dwarf to the 
nucleus of NGC 4438. 
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exhibit a X-ray counterpart. Two velocity components are 
observed in some pixels of the data-cube in the approach¬ 
ing side of the disk. One of the component has a redshifted 
kinematics which totally contrasts with the local, blue- 
shifted velocity field of the approaching half. Its velocity 
is however consistent with that found in nearby, diffuse, 
off-plane knots, suggesting that both could belong to a 
same morphological structure. The global kinematics of 
the extended filaments indicate that they are likely bound 
to the galaxy. 

(3) A tight relation between the presence of a small 
stellar tail aligned with a dust lane and an Ha filament 
South of the galaxy, and of the dwarf elliptical galaxy 
VCC 1040 is outlined. The gaseous filament suddenly 
changes its orientation at the tip of the stellar stream. 
The VCC 1040 redshift is comparable with velocities ob¬ 
served in the Ha filament. We suggest that a minor tidal 
interaction has occurred, or is perhaps occurring, between 
NGC 4438 and the dwarf galaxy, causing the disruption 
of the small galaxy and the formation of the stellar tail. 

Further 3D numerical simulations taking into account 
both tidal interactions and ram pressure stripping would 
be crucial to improve our understanding of the evolution 
of this complex system. 
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Fig. 7. (a). F450W WFPC2-H5T image showing VCC 1040 and the southern parts of NGC 4438. Full-line contours 
represent the Ha emission in NGC 4438, a dashed line contour an isophote of VCC 1040. Labels F2 to F4 refer to the Ha 
filaments. The opposite extremities of the stellar tail T1 are delineated by labels Tl. Numbers into brackets represent 
radial velocities in F4 and of VCC 1040 (in km s -1 ). A cross marks the location of the object 2MASXiJ1227422+125935 
(see text), (b). F450W-F814W colour index map. Dark shades represent dust lanes. The dark area to the East of 
VCC 1040 is not due to dust but to the sky background level. 











